Drude weight in systems with open boundary conditions 
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For finite systems, the real part of the conductivity is usually decomposed as the sum of a zero 
frequency delta peak and a finite frequency regular part. In studies with periodic boundary condi- 
tions, the Drude weight, i.e., the weight of the zero frequency delta peak, is found to be nonzero for 
integrable systems, even at very high temperatures, whereas it vanishes for generic (nonintegrable) 
systems. Paradoxically, for systems with open boundary conditions, it can be shown that the co- 
efficient of the zero frequency delta peak is identically zero for any finite system, regardless of its 
integrability. In order for the Drude weight to be a thermodynamically meaningful quantity, both 
kinds of boundary conditions should produce the same answer in the thermodynamic limit. We 
shed light on these issues by using analytical and numerical methods. 

PACS numbers: 05.60.Gg,72.10.Bg 
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Transport properties define materials as superconduc- 
tors, metals, or insulatorsi^ In one dimension, trans- 
port can help differentiate integrable from nonintegrable 
systems. Integrable systems, in general, display an infi- 
nite conductivity at all temperatures^ an effect that has 
been related to the strong influence of conservation laws, 
while at high temperatures, generic (nonintegrable) sys- 
tems are expected to exhibit a finite resistivity caused by 
umklapp scattering and inelastic collisions. 

Let us consider a finite system of length L (in units of 
the lattice constant); the real part of the conductivitj*iS 
can be written as2^£ 
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is the Drude weight (or charge stiffness) , the stress tensor 
operator T = —q e lim.fc_>o i[ J(fc), n(—k)], J is the current 
operator, J nm = (n\J\m) are its matrix elements, n(k) is 
the Fourier transform of the local density operator hj , q e 
the charge, and T = (T)/HL. The Boltzmann weight of 
a state |m) is denoted by p m = e~^ tm /Z , with e m being 
its energy, (3 = 1/fesT, and Z the partition function. 

In one dimension, and finite temperatures, the Drude 
weight can also be computed as^^ (assuming that there 
is no true superconductivity) 
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infinite volume object D^ = lini[ ; _ ! . 00 Dl = Mthl^oo D L , 
and a reg is the "regular" part of the conductivity. In sys- 
tems with periodic boundary conditions (pbe's), Dl and 
D* L are found to be nonzero in integrable systems at all 
temperatures, i.e., they display an infinite conductivity, 
whereas Dl and D* L are expected to exponentially vanish 
with the system size for nonintegrable systems4i£ 

A very different situation arises when one studies sys- 
tems with open boundary conditions (obe's). There, 
one can legitimately define* the position operator X — 
j fi>j , which satisfies the following commutation re- 
lations: J = i[X,H], where H is the Hamiltonian of 
the system and T = — i[X, J]. Since J nm = i(e m — 
e n ){n\X\m), one can see that D* L ([3]) is identically zero, 
and that the first and second terms in the right hand side 
of Eq. @ cancel each other, so that Dl is also identically 
zero. The vanishing of Dl and D* L for systems with obe's 
is true regardless of their intcgrabilityjii 

Considering the results above, one could question 
whether the Drude weight is a meaningful thermody- 
namic quantity, i.e., whether one would obtain the same 
result, measuring it in (quasi-) one-dimensional (ID) ex- 
periments involving closed rings or open chains. This 
is a valid concern since Z?oo is not a usual bulk ther- 
modynamic quantity. (For the latter, the equivalence of 
boundary conditions in the TL is obvious.) In an expan- 
sion of the energy of the system in powers of L, Doo at 
T = is proportional to a term that scales like 1/L~& 
Our goal in this work is to shed light on this puzzle. 

We consider ID systems of spinless fermions, 
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Equation ([I]) is true for all boundary conditions. In 
the thermodynamic limit (TL), it leads to the decom- 
position Re [(700(0-1)] = TTD oc S(huj) + a reg (u>), where the 



with nearest (t) and next-nearest (t') neighbor hopping, 
interacting with nearest (V) and next-nearest (V) neigh- 
bor repulsive potentials. The sum over j is appropriately 
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defined for pbc's and obc's. &■ (cj) are the creation (an- 
nihilation) operators for spinless fermions at a given site 
j, and hj = CjCj are the corresponding density operators. 
This model is known to be intcgrablc for t' = V' = 0. (It 
can be mapped onto the XXZspin 1/2 chain.) To under- 
stand the effects of obc's in more generic nonintegrable 
systems, we study systems with finite values of t' and V . 

The current and stress tensor operators for this model 
can be written as 



j = !|*^ a t 6 . +1 + H 
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As a starting point for our analysis, let us consider 
the case of noninteracting fermions. With pbc's, the 
single-particle eigenstates of the Hamiltonian arc plane 
l/VLJ2j e jfemJ 'cj|0) with energy e m = 
where k m = 2-Km/L and m = —L/2 + 



waves | m 
—2t cos(k 



l,...,L/2. In this case, \J nm \ 2 = [2q e t sin(fc m )] 2 J m ,„, 
i.e., the matrix elements of J are fully diagonal, and 
Dl = (T)/HL = q 2 (H)/hL. At finite temperatures, and 
for sufficiently large system sizes, D* L = Dl- 

For obc's, on the other hand, standing waves be- 
ing the noninteracting eigenstates, have the form 
|m) = y/2/(L + 1) J2j sin(fc m j)cj|0) with energy e TO = 
— 2t cos(fc m ), where k m = 7rm/(L + 1) and m = 1, . . . , L. 
The matrix elements of the current operator arc then 
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sin 2 (fc m ) sin 2 (fc„) 
[cos(/c m ) - cos(fc n )]' 



and vanish whenever (m — n) is an even number, i.e., only 
eigenstates with a different parity produce nonvanishing 
values of J nm . This expression shows that the matrix 
elements of J have a very interesting property: although 
\Jnm\ vanish for m = n, they attain their largest values 
for the smallest (odd) differences between m and n. For 
very large systems, the factor 1/L 2 produces vanishing 
values of J nm unless n ~ m. Rewriting k n = k m ± 
irl/(L +1), / = 1,3,... one can see that Re [cr(w)] is a 
sum of delta functions S[lo ± 2nt sm(k m )l / (L + 1)] whose 
weight decreases as ~ l/l 2 . 

In Fig. [2(a), we show Re [<t(o;)] (binned to obtain a 
smooth curve) for noninteracting particles in two sys- 
tems with obc's and different sizes. This figure shows 
that (i) broadened delta functions collapse to larger peaks 
(whose width reduces as the system size increases) situ- 
ated at io ~ 2nl/ L with Z = 1,3,..., and (ii) these peaks 
move toward smaller frequencies as L increases. From 
the previous analysis, one can conclude that although 
for finite open systems Dl and D* L are always zero, for 
L — * oo, a delta peak develops at lo = 0, but this time 
generated by the collapse of delta peaks that come from 
the so-called regular part of the conductivity. In addi- 
tion, from the sum rule Re [cr(w)] duj = Tr(T)/hL^- 
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— L=16 (pbc's) 

— L=18 (pbc's) 

- - L=16 (obc's) 

- ■ L=17 (obc's) 




FIG. 1: (color online). Regular part of Re [or (w)] for integrable 
systems (t' = V' = 0) in (a) the noninteracting case (V = 0), 
(b) V — 1.0, and [inset in (b)] the isotropic (Heisenberg) 
point V = 2.0. All of these results were obtained at a very 
high temperature f3 = 0.001 (infinite in the noninteracting 
case) and half-filling. The inset in (a) shows how A (see text) 
behaves as one approaches the TL. The values of V, j3, and 
ui are given in units of t. 



one obtains that the weight of such lo = peak is 
ttDoo = limi_ ) . 00 Tr(T)/HL = \imL^oo7rq 2 (H)/HL, which 
is identical to TrD^ as obtained from periodic systems 
since the energy is identical in both cases. 

At this point, one may wonder about the behavior of 
the finite frequency {lo ~ ±2irl/L, Z = 1,3,...) peaks as 
L is increased. It may happen that as L — > oo, (i) all the 
weight is concentrated in the lowest lo = ±1/L peaks, 
and the others disappear, or (ii) the weight is distributed 
among several peaks with different values of lo. To answer 
this question, we have studied the summed weight of the 
peaks at lo ~ ±2n/L as a function of increasing system 
size, a quantity we call A. Results are shown in the inset 
in Fig. Ufa). They confirm the scenario (ii) above since, 
as L — * oo, A saturates at around 80% of 7rr. Hence, 
the other peaks with I > 1 remain finite, and they are 
needed to account for the exact Drude weight in the TL. 

Moving away from the noninteracting case, but keeping 
the system integrable, we cannot make the corresponding 
analytical treatment, so we turn to a full exact diagonal- 
ization of finite chains. We perform calculations in the 
grand-canonical ensemble, and all results presented in 
what follows are obtained at half-filling (n = 0.5). 

For generic integrable systems with pbc's, the matrix 
elements of J are not diagonal anymore like in the non- 
interacting case. This means that weight moves away 
from the Drude peak and the regular part of Re [c(w)] 
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FIG. 2: (color online), (a) Scaling of the Drude weight (pbc's) 
and A (obc's) with the system size, for V = 1.0 and 2.0 and 
p = 0.001. (b) The Drude weight, computed as Dl and D* L , 
and F as a function of for the isotropic case (V = 2.0). 



becomes finite. This can be seen in Fig. (Ub) and its in- 
set, where we have plotted Re [er(w)] for two integrable 
systems with different values of V. The scaling of the 
Drude weight with system size, together with a "simple 
minded" linear extrapolation to the TL, are shown in 
Fig.UJa). As seen there, for both V = 1.0 and V = 2.0, 
we obtain a finite value of Doo £d£iL 

The results presented so far were obtained at a very 
high temperature ((3 = 0.001), where finite size effects arc 
the smallest, but transport properties are still nontrivial. 
In Fig. [21b), one can see that for T > 10, the actual 
value of T is not essential since quantities such as Dl,T, 
D L T, and TT become almost independent of T. Figure 
[2Kb) also shows that, for the considered temperatures and 
system sizes, V is independent of the system size, while 
Dl and D* L still exhibit finite size effects. It is important 
to notice that albeit Dl and D* L are identical for any 
given system size at high temperatures, finite size effects 
build differences between these two ways of computing 
the Drude weight at lower temperatures £ 

After reviewing the pbc case, we can now analyze the 
effects of obc's on more generic integrable systems. In 
Fig. |Tfb) and its inset, we have plotted Re for 
chains with obc's, and two different values of V, together 
with the results for pbc's. One can clearly see there that 
large peaks develop in the obc data over the pbc results, 
and that these peaks move toward lower frequencies as 
L is increased 4 Actually, for V = 1.0, one can see that 
the two largest obc peaks are located at u> ~ 2irl/L, 




FIG. 3: (color online). Regular part of Re [<r(tj)] for nonintc- 
grable systems with pbc's (L = 17) and (inset) obc's (L = 16) 
as one moves away from the isotropic integrable point by in- 
creasing the value of t . In all cases, V — 2.0 and {3 = 0.001. 



I = 1 and 3, similar to the noninteracting case. These re- 
sults strongly suggest that in the TL, the finite frequency 
peaks for obc's will collapse into a single ui — Drude 
peak like the one obtained for pbc's. 

In Fig. [D^a) , we show how A behaves with increasing 
system size for V = 1.0 and V = 2.0. [A is computed in 
this case as two times the area over the dotted lines in 
Fig. [Ub).] In both cases, the behavior of A is consistent 
with extrapolated u> = peaks with a finite weight in 
the TL. Like for the noninteracting case, extrapolating 
A to L — * oo does not reproduce the value of Z?oc ob- 
tained for systems with pbc's. This, we infer, is due to 
the weight distributed over peaks with higher frequen- 
cies, which, following the noninteracting case, should all 
collapse to co = as L — * oo. The extrapolations in Fig. 

show that the relative weight of peaks with higher 
cu in chains with obc's increases as one departs from the 
noninteracting case. 

Having shed light on integrable systems, we now turn 
to the nonintegrable case. As mentioned before, in sys- 
tems with pbc's, breaking intcgrability is expected to 
produce an exponentially (with the system size) vanish- 
ing Drude weight. We first consider (Fig. [3]) the case 
in which in Eq. Q integrability is broken by t' , which 
also breaks the particle-hole symmetry of the integrable 
model. In Fig. [3l one can see that, even for small finite 
systems, introducing t' has dramatic effects in Re [<r(u;)]. 
A peak develops at low frequencies. Given the sum rule 
for Re [o"(w)] , such a peak can be related to the disappear- 
ance of the Drude peak and the transfer of its weight to 
finite frequencies. Hence, the system acquires a finite dc 
conductivity, which decreases with increasing t' (Fig. [3]). 

For small finite systems with obc's, on the other hand, 
nothing dramatic should happen when a small t' is in- 
troduced. This is because in the integrable case, there 
is no zero frequency delta peak, but, instead, finite fre- 
quency peaks that are already present in the regular part 
of Re [cr(w)] . As seen in the inset in Fig.[31 adding t' to in- 
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FIG. 4: (color online). Regular part of Re [<r(w)] for nonin- 
tegrable chains with different system sizes, (a) t-t'-V model 
with t' = 0.5, (b) t-V-V model with V' = 1.0, and [inset 
in (b)] t-t'-V model for a system of hard-core bosons with 
t' = 0.5. The inset in (a) shows the low-frequency part of 
Re [cr(o;)] for the fermionic t-t'-V model together with a non- 
linear fit (see text). In all cases, V = 2.0, /3 = 0.001, and t' , 
V, V , /3, and ui are given in units of t. 

tcgrablc chains with obc's only reduces the height of the 
lowest frequency peak in Re [c(w)] and, with increasing 
i', its weight moves toward higher frequencies. 

An apparent difference between Re in systems 

with pbc's and systems with obc's (Fig. [3]) is that, while 
the former exhibit a finite dc conductivity, the latter 
exhibit a vanishing one. The vanishing of the dc con- 
ductivity for finite systems with obc's is understandable 
by using an analogy with pbc's with a small but finite 
momentum q. For q =/= 0, we know from the conserva- 
tion law [wh{q,uj) ~ qJ(q,u>)] that the current must 



vanish in the dc limits One expects, in this case, to 
see diffusion, i.e., a(q 7 uj) = ct(0,cj)/[1 — iDq 2 /ui] with 
q ~ tt/L for obc's. This means that at very low fre- 
quencies for obc's, Re [cr(w)] ~ uj 2 . As shown in the 
inset in Fig. Ufa), our data are crudely consistent with 
that. There, we also show the result of a nonlinear fit to 
Rc [cr(u)] — J2a=i 6 CaW a , where we find that C\ ~ 0. 

From the previous analysis, one expects the region of 
w's over which Re [cr(w)] — > to reduce with increasing L 
and, eventually, in the TL, to recover the k = result ob- 
tained for pbc's. In the main panel of Fig.fJJa), we depict 
results for pbc's and obc's in systems with two different 
sizes. There, one can see that, indeed, for obc's, the low- 
frequency region with decreasing conductivity moves to 
lower frequencies with increasing system size. In the TL, 
the usual belief is that long time tails in the autocorrela- 
tion of the current will ultimately take over, asymptoti- 
cally leading to Re ~ A — By/\uj\£ 

We should stress that our conclusions above are valid 
for generic nonintcgrablc systems, i.e., they are not lim- 
ited to the t-t'-V model presented in Figs. [3] and[4|a). 
For example, in Fig. 0Jb) we show that similar results 
are obtained when t'=0 but V'^0 (t-V-V model, main 
panel), and V'=0, t'^0 (t-t'-V model), but for a system 
of hard-core bosons (inset) . The last two models preserve 
the particle-hole symmetry present in the intcgrable case, 
so our conclusions are also independent of its presence or 
absence in nonintegrable systems. 

In summary, we have argued that even though the real 
part of the conductivity in finite systems with obc's is 
qualitatively different from that of systems with pbc's, 
they both have a common thermodynamic limit. For 
integrable systems, this means that there is an infinite 
dc conductivity, characterized by a finite Drude weight. 
On the other hand, for nonintegrable systems, the dc 
conductivity is finite and, in our ID systems, it decreases 
as one moves away from the integrable point . 
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